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0^ ' The existence of edge states and zero energy modes in vortex cores is a hallmark of topologically 

' nontrivial phases realized in various condensed matter systems such as the fractional quantum Hall 

, states, p + ip superconductors, and Z2 insulators (quantum spin Hall state). We examine this 

scenario for two dimensional noncentrosymmetric superconductors which allow the parity-mixing 
pH ' of Cooper pairs. It is found that even when the s-wave pairing gap is nonzero, provided that 

, the superconducting gap of spin-triplet pairs is larger than that of spin-singlet pairs, gapless edge 

states and zero energy Majorana modes in vortex cores emerge, characterizing topological order. 
Furthermore, it is shown that for Rashba superconductors, the quantum spin Hall effect produced by 
gapless edge states exists even under an applied magnetic field which breaks time-reversal symmetry 
provided that the field direction is perpendicular to the propagating direction of the edge modes. 
This result making a sharp contrast to the Z2 insulator is due to an accidental symmetry inherent 
H ' in the Rashba model. It is also demonstrated that in the case with magnetic fields, the non-Abelian 

^ , statistics of vortices is possible under a particular but realistic condition. 
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I. INTRODUCTION 



Topological phases of matter are quantum many-body states with topologically non-trivial structures of the Hilbert 
spaces [1]. They are characterized by the existence of both topologically-protected gapless modes on the edges and 
' ^ I a bulk energy gap which separates the ground state and excited states [2-6]. Recently, topological phases are of 
^ i interest in connection with wide ranges of subjects in condensed matter physics, such as the quantum Hall effect [7], 
Q _ p -\- ip superconductors [3-5, 8-10], and Z2 topological insulators (quantum spin Hall effect) [11-13]. The gapless 
edge states are topologically stable against the perturbations that do not break symmetries of systems, and play a 
crucial role in transport properties such as quantum (spin) Hall effects, which is, currently, inspiring applications to 
^ spintronics [2, 13, 14]. The total number of topologically-protected edge modes in a given system is associated with 
topological numbers, i.e. the TKNN number (the first Chern number of the U(l) bundle describing the many-body 
\^ wave function) for systems without time-reversal symmetry [6, 7], and the Z2 invariant (the parity of the spin-resolved 
00 ■ TKNN number) for time-reversal invariant systems [11]. 

The presence of low-energy edge modes also implies the fractionalization of quasiparticles in the bulk [15]. For 
instance, in a vortex core of a spinless p -\- ip superconductor, there is a zero energy mode which is described by a 
Majorana fermion, i.e. a half of a conventional fermion [16]. A vortex with a Majorana fermion is a quasiparticle 
00 ' obeying the non-Abelian statistics [8-10, 17, 18]: The braiding of vortices with Majorana fermions gives rise to the 
■ superposition of the degenerate many-body ground states. Owing to this entangled character, non-Abelian anyons can 
be utilized for the construction of fault-tolerant quantum computers [18-22]. When there are odd number of vortices 
with Majorana fermions, the edge state must be also Majorana, since an isolated Majorana fermion is unphysical, 
and should be paired with another Majorana fermion [23]. In this sense, the Majorana edge state is a concomitant of 
a zero energy Majorana mode in a vortex core, and the existence of both of them characterizes the topological order. 

A more precise argument is given in terms of the low-energy effective theory. The low-energy effective theory for 
topological phases is the Chern-Simons theory, and both edge states and topological quasiparticles in the bulk are 
described by the corresponding conformal field theory. In the case with the level-/c SU(2) symmetry, it is the SU(2)fc 
Wess-Zumino-Witten theory, which is decomposed into the U(l) gaussian theory with the central charge c = 1 and the 
Zk parafermion theory with c = "^^^^2"^ [3~5]. For superconductors, where quasiparticles are superpositions of particles 
and holes, the U(l) gaussian part describing fractional charges for the fractional quantum Hall state is absent, and 
thus, the low-energy effective theory for edge states and quasiparticles is the Zk parafermion theory. The topological 
phase of p + ip superconductors corresponds to fc = 2, i.e. the Ising conformal field theory, the operator content of 
which includes Majorana fermions [18, 24]. 

The search for possible realization of topological phases is an intriguing and challenging issue involving the devel- 
opment of novel concepts for quantum condensed states as well as potential technological applications. In this paper. 
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we demonstrate that two dimensional (2D) noncentrosymmetric superconductors (NCS) under certain circumstances 
provide another candidate for the physical realization of a topological phase, which supports the existence of edge 
states and zero energy modes in vortex cores. In NCS, asymmetric spin-orbit (SO) interactions play crucial roles in 
various exotic superconducting phenomena such as parity-mixing of Cooper pairs, magnetoelectric effects, and helical 
vortex phases [25-33]. It is found that the asymmetric SO interaction is also a key to realize topological phases in 
NCS. We investigate the edge states and the zero energy vortex core states in NCS by using both numerical and 
analytical methods, and verify the condition for the realization of topological phases in NCS. 

The main results of this paper are as follows. For 2D NCS with the admixture of s-wave pairing and p-wave pairing, 
as long as the p-wavc gap is larger than the s-wave gap, topological phases are realized. The classification of the 
topological phases in NCS is completed by examining topological numbers. In the absence of magnetic fields, a 
topological phase with two gapless edge modes emerges, which leads to the quantum spin Hall effect; i.e. when an 
electric field or a temperature gradient is applied, a spin Hall current carried by the edge state flows in the direction 
perpendicular to the applied external field. Moreover, it is shown that in the case of the Rashba SO interaction, the 
gapless edge modes are stable against weak magnetic field applied perpendicular to the direction in which the edge 
modes propagate. This is a bit surprising because the magnetic field which breaks time-reversal symmetry flaws the 
Zi classiflcation of the topological phase, and also the TKNN number is still zero for such a weak magnetic field. 
In fact, the stability of this topological phase is due to an accidental symmetry for particular symmetry points in 
the Brillouin Zone specific to the Rashba model, which is characterized by another topological number, "winding 
number" . We will clarify the condition for the nonzero winding number in this paper. In addition to these spin Hall 
states, there is also a topological phase with the nonzero TKNN number, which is analogous to the quantum Hall 
state, and realized for a particular electron density with a magnetic field. The implication of these gapless edge states 
for experimental observations will be also discussed. 

Furthermore, we examine the vortex core state for these topological phases by solving the Bogoliubov-de Gennes 
equations. On the basis of both numerical and analytical methods, it is found that zero energy Majorana fermion 
modes exist in the vortex cores. It is also verified how the non-Abelian statistics of vortices is realized in NCS. The 
non-Abelian statistics of vortices in superconductors has been considered so far for the chiral -p^ state. [8, 9] In 
the case of spinless p + ip state, there is only one zero energy Majorana mode in a vortex, which is a non-Abelian 
anyon. In the case of spinful p + ip state, there are two Majorana fermions in a vortex corresponding to spin up and 
spin down states, which form a usual complex fermion. The non-Abelian statistics is not realized for this situation. 
As Ivanov elucidated, when there is a half quantum vortex which suppresses one of two Majorana fermions in a 
vortex core, vortices behave as non-Abelian anyons [9]. For NCS, however, this scenario is not applicable, because the 
half-quantum vortex is a texture of the d-vector of p-wave pairing, and the d-vector in NCS is locked by the strong 
SO interaction. As proposed by one of the authors [34], when one tunes a chemical potential so as that the Fermi level 
crosses the F point in the Brillouin zone at which the electron band is a Dirac cone, and applies a Zeeman field parallel 
to the z-direction, a Majorana fermion associated with the Dirac cone is eliminated, and only one Majorana fermion 
survives in a vortex core, which makes the non-Abelian statistics of vortices possible. Wo examine this strategy by 
using an analytical approach based on the index theorem for zero modes developed by Tewari, Sarma and Lee in the 
case of p + ip superconductors [35]. It is found that when the spin-triplet gap is larger than the spin-singlet gap, the 
non-Abelian statistics of vortices is realized for the situation mentioned above. 

It would be useful to comment on some recent other studies related to the present paper. 1) Some topological 
properties of the Rashba type NCS, such as the nodal structure of the gap function and the edge state, were discussed 
in [36, 37]. While these literatures assumed the time-reversal symmetry and were based on the Z2 topological number, 
we complete the topological classification by taking into account the TKNN number simultaneously. We also find that 
the "winding number" mentioned above is useful to characterize topological phases of the Rashba superconductors. 2) 
In [34] , a possibility of zero energy states in vortex cores for purely s-wave noncentrosymmetric superconductors under 
strong magnetic fields was suggested. Later, unfortunately, it was turned out that the zero energy vortex core state 
found in [34] is spurious, and there is no zero energy state in purely s-wave cases. The misleading conclusion in [34] 
is due to the erroneous assumption that the unitary transformation which diagonalizes the spin index is commutative 
with the center of mass coordinate of Cooper pairs raised by the presence of vortices. This assumption is valid only 
when the quasiclassical approximation is justified. This is not the case for the issue of zero energy vortex core states. 
This point will be discussed in more detail in Sec. IV B of the current paper. 3) Independently, Lu and Yip elucidated 
that in NCS with parity-mixing, the zero energy vortex core state is possible only when the gap of the spin-triplet 
pairs is larger than that of the spin-singlet pairs [38]. They examined the condition for zero energy modes of the 
Bogoliubov-de Gennes (BdG) equations without obtaining explicit solutions. However, to investigate the possible 
realization of the non-Abelian statistics of vortices, it is desirable to derive the explicit zero energy solutions which 
enable us to see whether a zero energy mode in a vortex core is a Majorana fermion or a usual complex fermion. This 
issue is addressed in the current paper. 

The organization of this paper is as follows. In Sec. II, we start with a general topological argument, which enables 
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us to understand edge and vortex core states clearly, providing the complete classification of topological phases in 
NCS. In Sec. Ill, gapless edge states in NCS are investigated by using the general topological argument and numerical 
methods. The implication for transport properties associated with edge states is also discussed. In Sec. IV, we explore 
zero energy vortex core states for NCS on the basis of numerical solutions for the Bogoliubov-de Gennes equations and 
an analytical approach based on the index theorem for zero energy states, and clarify the condition for the realization 
of the non-Abelian statistics of vortices. We conclude in Sec. V with a summary of our results. 

II. TOPOLOGICAL NUMBERS FOR NONCENTROSYMMETRIC SUPERCONDUCTORS 

In this paper, we consider type II NCS with Rashba-type SO interaction [39] in two dimensions. For concreteness, 
we define our model on the square lattice, though the following argument does not rely on the particular choice of 
the crystal structure. Then the model Hamiltonian is 

+ \ Y. ^-'(fc)cLcLfc,. + ^ ^ A;,,(fc)c_fe,Cfc,-, (1) 

kcFCF' ktrtr' 

where c|,^ (cka) is a creation (an annihilation) operator for an electron with momentum k = {kx,ky), spin a. The 
energy band dispersion is Sk = —2t{cosk,j. + cosfcj,) — /i with the hopping parameter t and the chemical potential /i, 
and the Rashba SO coupling is Co{k) = {smky, — sinkx)' For simplicity, we assume that t > and < 0. Because of 
parity mixing of Cooper pairs, the gap function A(fc) has both a spin-triplet component d{k) and a spin-singlet one 
'(/'(fe) at the same time, A(A;) = it/j{k)ay + id{k)a(jy. Due to the strong SO coupling, the spin-triplet component d{k) 
is aligned with the Rashba coupling, d{k) = AtCoik) [29]. For the spin-singlet component ip, we assume a s-wave 
pairing, il){k) = Ag. The amplitudes At,s are chosen as real. The Zeeman coupling Mb-^z X^fc(c][.|Cfc| — cl.^Cki) with 
Hz a magnetic field in the z direction has been also introduced for later use. 

Before going to study topological properties of the system, we first examine the bulk spectrum of the system. 
Topological nature of the system changes only when the gap of the bulk spectrum closes. The bulk spectrum E{k) of 
the system is obtained by diagonalizing the following matrix, 

TT(i,\ - ( '^k- MB-H"zcr2 + aCo{k) ■ a iA^ay + iAtCo{k) ■ acjy \ , , 

^^•^l-y -iA,ay - iAtCo{k)ay a -£fc + piY^E^a, + a£o(fe) ■ a* ) ' 

and we have 



E{k) = ±]jel + (a2 + A2)£o(fe)2 + f^^H^^ + ± 2^ {ska + A^A.fCoikf + (4 + A2)^|if|. (3) 
The gap of the system closes only when 

4 + {a' + ADCoikf + nlHl + = 2sJ{eka + A,A,f Ca{kf + {el + A2)^|if|, (4) 
which is equivalent to 

el + Al=iilHl + {a'' + Al)Co{kf, ekA,Co{k) = A,aCo{k). (5) 
When At ^ 0, (5) is met either when 

= t;^' (l + ^) (A? A(fc)^ - A^) + t^lHl = 0. (6) 



or 

-2 I a2 ..2 tj2 



6i + Ai = ij,iHt, £o(fe) = 0. (7) 
In the absence of the magnetic field, only Eqs.(6) can be met and they are rewritten in simpler forms, 

el = a'Co{kf, AlC,{kf=Al. (8) 
Topological nature of the system does not change unless (6) or (7) (or (8) when Hz = 0) is satisfied. 
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A. Z2 topological number 



When Hz = 0, the system is time-reversal invariant, and the topological property is characterized by the Z2 
invariant [11, 40-44], [59]. We will show that if the spin- triplet pairing is stronger than the spin-singlet one, the Z2 
number is non-trivial. 

To calculate the Z2 number, we adiabatically deform the Hamiltonian of the system without gap closing. This 

process does not change the Z2 topological number, since it changes only when the gap closes. From (8), it is 
found that if the spin-triplet amplitude AtjCo(fc) is larger than the spin-singlet one Ag on the Fermi surface given by 
Sk — oiCo{k), we can take Ag 0, then a without gap closing. (If ej, = at one of the time- reversal momenta 
fe = (0, 0), (tt, 0), (0, tt), (tt, tt), the gap closes when Ag = 0. However, this undesirable gap closing can be avoided by 
changing ^ or t slightly.) Thus, its Z2 number is the same as that of the pure spin-triplet SC with d{k) = AtCo{k). 

As was shown in [45], the Z2 topological number (— 1)"^ for a pure spin-triplet SC is determined by its dispersion 
£fe in the normal state, 

(-1)-= n sg^^fe- (9) 

fe=(0,0),(0,7r),(7r,0),(7r,7r) 

From this formula, it is easily shown that the Z2 number is always non-zero (mod 2) for the square lattice system. 
Therefore, from the argument above, if the spin-triplet pairing is stronger than the singlet one, the NCS is a topological 
insulator with a non-trivial Z2 number. This topological superconductor belongs to the same class as one discussed 
in [37, 46, 47]. 



B. TKNN number 



In the presence of a magnetic ficild, the timc-rcvcrsal invariance is broken. The Z2 topological number is no longer 
meaningful, and the TKNN number plays a central role in topological nature of the system instead. In the following 
argument, we consider only the Zeeman effect of the magnetic field, neglecting the orbital depairing eff'ect. When the 
energy scale of the SO interaction is sufficiently larger than the Zeeman energy scale, the Pauli depairing effect docs 
not exit for magnetic fields parallel to the z-axis, and the structure of the d-vector for the spin-triplet pairs is not 
affected by the Zeeman effect [29, 33]. We assume such situations. The effect of the orbital depairing effect will be 
discussed at the end of this subsection. 

In order to obtain a non-zero TKNN number, the magnetic field should be large enough to have a gap closing. 
Otherwise the TKNN number must be zero since the system is adiabatically connected to the case with = 0. 
(When Hz = 0, the TKNN number becomes zero because of time-reversal invariance.) When the spin-triplet pairing 
is stronger than the singlet one, (6) is never met in the presence of the magnetic field, and a gap closing occurs only 
when (7) is satisfied. Since Co{k) is zero at fe = (0,0), (0, tt), (tt, 0), (tt, tt), it is found that we have a gap closing if 
one of the following equations is satisfied, 

{-4t-nf + A',=nlHl /x2 + A2=^|Ff, {4t-nf + Al=nlHl (10) 

From this, it is found that we have a two classes of phase diagram illustrated in Fig.l. a) For —At < n < —2t, the 
first phase transition occurs at hbHz = \/ (4i -I- /i)^ + A^, and the second one at hbHz — ^ tjfi + A^, and the third 
one at ijlbHz = \/ (4f — \xf- -|- A^. b) For —It < /x < 0, the first phase transition occurs at \x^Hz = ■\//i^"+Af, and 
the second one at [i^Hz = y^(4F+7tp"-i-Af, and the third one at \x^Hz = 1/ (4t — /i)^ -|- A^. The patterns of the 
TKNN numbers differ for these parameter regions. 

The TKNN number for each phase is evaluated by using the adiabatic deformation of the Hamiltonian in a similar 
manner as Sec. II A. As is evident from Fig.l, we can take Ag = a = without crossing the phase boundary. For 
As = a = 0, the up-spin electron and down-spin electron are decoupled from each other, and the Hamiltonian reduces 
to a pair of 2 x 2 matrices. 



jj (k) = ( ^k-mHz -iAi{smkx - isinky)\ 
' yiAt{smkx+isva.ky) -£k + I^BHz J 



i At (sin fca; - i sin /cy ) -e^- [i^H, 
Using these 2x2 matrices, we can evaluate the TKNN integers. By rewriting H„„{k) (cr =t, i) as H„„{k) = R„{k)-a 
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FIG. 1: Schematic picture of phase diagram of 2D NCS with Rashba couphng. 



the TKNN number /tknn is given by 



-^TKNN 




(12) 



with R„ = R„/\R„\. 

In Table. I, we summarize the obtained TKNN numbers. The phases IV and IV' are trivial band insulators, and 
of no interest. The other phases are in topologically nontrivial states. For the phases II, III, and III', there is only 
one band which crosses the Fermi level, and is associated with the nonzero TKNN number. The sign of the TKNN 
number depends on the curvature of the Fermi surface as well as the chirality of the gap function, since it is related 
to the Hall conductivity. The gap functions for the phases II, III, and III' possess the same chirality, and hence, for 
the positive curvature of the Fermi surface, /tknn > (the phase II), and for the negative curvature, /tknn < 
(the phases III and IIP). In the phase IF, there are two bands with opposite signs of the curvatures of the Fermi 
surfaces. Also, the gap functions for these bands have opposite chiralities. Thus, each of the two bands contributes 
to the TKNN number equal to -1, leading to the total TKNN number /tknn = -2. The phases II, III, IF, and III' 
are analogous to the quantum Hall state. 

For the phases I and I', although both the TKNN number and the Z2 index are zero, there exist topological orders 
with two gapless edge modes and Majorana fermions in vortex cores, as will be shown in Sec. HI and Sec. IV, and 
their topological nature is similar to that of the Z2 insulator. In the next subsection, it will be clarified that these 
phases are characterized by another topological number, i.e. winding number. 

In the argument above, the orbital depairing effect of magnetic fields is ignored. For typical superconductors, the 
orbital limiting field is H°^^ ^ Ep{A/ Ep)'^ / fiB ^ t/f^B with the Fermi energy Ep. Thus, superconductivity does not 
survive in the strong magnetic field regions. III and HI', though the Pauli depairing effect is completely suppressed 
by the strong asymmetric SO interaction as mentioned before. (To avoid confusion, we would like to stress again that 
the topological phases obtained above by putting a — can be deformed into topologically equivalent states with 
nonzero large a, unless the bulk gap closes.) On the other hand, the realization of the phases I, II, P, and IP in the 
weak field regions is feasible. 



As was mentioned in the previous subsection, the NCS with the Rashba coupling has an additional topological 
number, in addition to the Z2 and TKNN numbers. Unlike the usual topological number, the additional topological 
number is defined only for particular values of momentum, but it is also useful to understand properties of edge states 
and vortex core states, which are specific to the 2D NCS, in the presence of a magnetic field. 

To define the topological number, let us consider the particle-hole symmetry of the Hamiltonian, 



C. Winding number 




(13) 



For ky = or ky = tt, we have H{—k)* = H{k), thus the Hamiltonian anti-commutes with P, {T,H{k)} = 0. This 
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TABLE I; The TKNN integer Jtknn and the winding number I{ky) for 2D NCS with Rashba coupling, 
implies that if we take the basis where T is diagonal F = diag(l, 1,-1,-1), then the Hamiltonian becomes off-diagonal 



q{k) 



gt(fe) 



H{k) = 

for these values of ky. Using q{k), we can define the topological number as [48], 

I{ky) = 4^ dkMQ~\k)dkMk) - q^-\k)dkj{k)), {ky 
An explicit calculation shows that q{k) is given by 

q{k)-- 



0, tt). 



£fe - /iB-ffz +iAtsinfca: -Ag + iasinfca: 
As - ia sin k^ Sk + Ms-f^z + i^t sin k^ 



(14) 



(15) 



(16) 



From this, wc can calculate I{ky) {ky = 0, tt) for each phase in FigT. The result is summarized in Table J. 

The additional topological number I{ky) is accidental and very sensitive to the direction of the magnetic field. 
While I{ky) remains well-defined even in the presence of a magnetic field in the a;-direction, it becomes meaningless if 
wc apply a magnetic field in the y-direction, Hy^ since Hy breaks the relation H{—k)* = H{k) for ky = 0, tt. From the 
bulk-edge correspondence, this means that edge states for Rashba type NCS are also very sensitive to the direction 
of the magnetic field, which will be confirmed numerically in the next section. 

III. EDGE STATES IN NONCENTROSYMMETRIC SUPERCONDUCTORS 

In this section, we investigate edge states for the 2D NCS numerically. Prom the bulk edge correspondence, a 
non-trivial bulk topological number- implies gapless edge states. This will be confirmed in this section. Experimental 
implications for the gapless edge states are also discussed in Sec. Ill C. 



A. Without a magnetic field 

Let us first study edge states for the 2D NCS in the absence of magnetic field. Consider the lattice version of the 
Hamiltonian (1) 



(i,J>,<T,cr' 



i 



(17) 
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FIG. 2: The energy spectra of the 2D NCS with edges &t = and = 50 in the absence of magnetic field. Here ky denotes 
the momentum in the j/-direction. We take t = 1, jj, = —3, a = 0.6. a) NCS with dominating p-wave paring. At = 0.6 and 
As = 0.1. b) NCS with purely s-wave paring. At = and Ag = 0.6. 

where i = {ix.iy) denotes a site on the square lattice, rij a unit vector from a site i to a site j. The sum J2{ij) 
is taken between the nearest neighbor sites. In this subsection, we suppose = 0. Consider the system with two 
edges at = and = , and put the periodic boundary condition in the y direction. By solving numerically the 
energy spectrum as a func;tion of the momentum ky in the y direction, edge states for NCS are studied. 

As was shown in the previous section, the Z2 topological number is non-zero if the gap of the spin-triplet pairs are 
larger than that of the singlet. Therefore, from the bulk-edge correspondence, there should always exist gapless edges 
if the spin-triplet pairs dominates the superconductivity. In Fig. 2 a), we show the energy spectrum of the 2D NCS 
with edges. It is found that there exist gapless edge states in the bulk gap. The gapless edges states form a Kramers 
pair. 

For comparison, we also illustrate the energy spectrum for the 2D NCS with purely ,s-wavc paring in Fig. 2 b). As is 
seen clearly, no edge state is obtained. This is also consistent with the trivial Z2 number of the purely s-wave paring. 



Let us now examine edge states in the case with a magnetic field. As is shown in the previous section, there exists 
a variety of topological phases characterized by the topological numbers. 

In Figs. 3 and 4, we illustrate the energy spectra of 2D NCS with edges for various topological phases. All phases 
have a bulk gap, and some of them have gapless edge states corresponding to the non-trivial topological numbers 
in Table I. It is found that for ky with nonzero I{ky) {ky = 0,7r), a zero energy edge state always appears, and 
the number of zero energy edge states coincides with the absolute value of I{ky). We also find that a phase with a 
non-zero /tknn has a edge state with the total chirality /tknn ■ These results are also consistent with the bulk-edge 
correspondence. 

Wc also notice that the gapless edge states in the phases I and I' are very sensitive to the direction of the applied 
magnetic field. As seen in Fig. 5, while the gapless edge states are stable imder a magnetic field in the x-direction, 
they become unstable under a small magnetic field in the y-direction. This bc;havior is naturally understood by the 
sensitivity of the definition of I{ky) to the direction of the magnetic field, which was mentioned in the previous section: 
For the phases I and F, the gapless edge states are ensured by I{ky), but in the case with non-zero Hy its existence is 
no longer protected since the winding number becomes ill-defined. As a result, the magnetic field Hy along the edge 
causes a tiny gap of the order 0{^BHy) for the edge states. 



B. With a magnetic field 



C. Transport phenomena associated with edge states 



The existence of gapless edge states revealed by the previous subsections implies that transport phenomena asso- 
ciated with the edge states are possible in analogy with the quantum Hall state and the Z2 topological insulator. 
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HbH, = 1.5, III) hbH, = 5, and IV) hbH^ = 8. 



Here, we discuss such phenomena in NCS. Transport properties inherent in edge states can be probed experimentally 
by using the Hall bar geometry considered before for the detection of edge states of the quantum (spin) Hall effect 
depicted in Fig. 6 [49-51]. Since our systems are superconductors, it is important for the experimental detection 
to discriminate between contributions from supercurrents and currents carried by edge states. A simple approach 
suitable for this purpose is to use thermal transport measurements. To suppress contributions from the Bogoliubov 
quasiparticles in the bulk, we assume that temperature is sufhciently lower than the superconducting gap, T <C A, 
and also the superconducting gap does not have nodes where the gap vanishes. The thermal conductance for heat 
currents is defined by = /^/(AT)i4 where is a thermal current between contacts i and j in Fig. 6, and 
(AT)y is the temperature difference between these contacts. In contrast to the conductance for electric currents, the 
thermal conductance is not quantized but depends on temperature T. The T-dependence of governed by edge 
states obeys a power law oc T, which can be distinguished from contributions from the bulk quasiparticles which decay 
like ~ exp(— A/T). Furthermore, as discussed in [51] in the case of the quantum spin Hall effect, in a six terminal 
measurement, there is no temperature difference between contacts 2 and 3 (or 5 and 6), because the edge current is 
dissipationless. 

A direct probe of spin Hall current carried by edge states may be also possible by measuring magnetization due to 

spin accumulation at contacts as discussed for the spin Hall effect [14]. Bulk supercurrents carried by Cooper pairs 
do not contribute to spin Hall currents even for a spin-triplet pairing state, and hence the spin Hall current between 
contacts 3 and 5 induced by a longitudinal voltage or temperature difference between contacts 1 and 4 is governed by 
edge states. 

A more remarkable effect due to edge states is the existence of the non-local transport [49, 50]. The non-local 
conductance is given by a heat current between contacts 3 and 5 divided by the temperature difference between 2 



9 




FIG. 5: The dependence of gapless edge states in the phase I on the direction of the magnetic field. We take t = 1, fj, = —3, 
a = 0.6, At — 0.6 and As = 0, 1 a) fiB-Hx = 0.15, usHy — and fisHz = 0. b) fisHx = 0, usHy = 0.15 and hbHz = 0. A 
similar dependence is obtained for gapless edge states in the phase I'. 
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FIG. 6: Hall bar with six terminals used for the detection of edge states. Solid and dotted gray lines with arrows represents 
propagating edge modes with different spin chirality corresponding to two SO split bands. These modes propagate in opposite 
directions. 

and 6 in Fig. 6, GJ^l = l'iz/{^T)2Q- If contacts 3 and 5 are well separated from 2 and 6, the nonzero G^^^ can not 
be explained by the bulk quasiparticles, which provides a direct evidence for the existence of current-carrying edge 
states. 

In the case with a sufficiently weak magnetic field, the spin Hall effect still exists as long as the direction of the 
magnetic field is perpendicular to the propagating direction of the edge states, because of the accidental symmetry 
of the Rashba model as discussed in Sees. 11 C and IIIB. If the magnetic field is tilted, and the field component 
along the propagating direction is nonzero, a gap opens in the energy spectrum of the edge states; this leads to the 
suppression of the spin Hall current which can be observed as a drop in the temperature difference between contacts. 

The gapless edge states are also observed experimentally as a zero bias peak of tunneling current [52] . Our results 
suggest that structures of the zero bias peak is very sensitive to the direction of the magnetic field. Under a magnetic 
field perpendicular to the edge direction, the zero bias peak is observed, but under a magnetic field along the edge, a 
splitting of the conduction peak is observed due to a tiny gap of the gapless edge states. 

The transport properties considered above are experimentally observable for 2D Rashba NCS with At > Ag which 
do not possess gap-nodes. 

IV. MAJORANA ZERO ENERGY MODES IN VORTEX CORES AND THE NON-ABELIAN 

STATISTICS 

In this section, we explore zero energy states in vortex cores for 2D NCS on the basis of both numerical and 
analytical approaches. We study the case with purely p-wave paring and the case with the admixture of s-wave and 
p-wave parings, respectively. It is assumed that no node exists in the superconducting gap because the existence of a 
full energy gap in the bulk is crucial for the stability of topological phases. As mentioned in the previous sections, as 
long as the gap of p-wavc pairing is larger than that of s-wave pairing, the state is topologically equivalent to the purely 
p-wavc pairing state, exhibiting topological nontriviality. In the following, we confirm this topological argument by 
obtaining explicitly the zero energy solutions of vortex cores. The condition for the non-Abelian statistics of vortices 
is also clarified on the basis of the explicit solutions. 

A. Numerical analysis of BdG equations 

In this subsection, we analyze the low energy states in a vortex core of NCS by using numerical methods. For 
this purpose, we consider the following two dimensional tight-binding model for a Rashba superconductor with the 
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FIG. 7: The density of quasiparticles for the purely p wave state plotted on the 81 x 81 xy-plane for some low energies E. 
The top and bottom panels are, respectively, the plots of |u-f(r)p + [u|(r)p and |u|(r)p + |f|(r')p for a = t, At = 0.05t, 
/XB^^z = 0, A* = -3.75t. E = 4.87 x 10""'* (left), E = 7.67 x 10"''* (middle), E = 1.023 x 10"^* (right). For E = 4.87 x 10""* 
and E = 1.023 x lO^'^t, edge states appear. 

admixture of s-wave and p-wave pairings, 
+A.^(e*cl^4,+h.c.) 

i 

-1a, ^(e^^-f ct,4,, - ct,4,, + e^^-f + ^e'^^i + h.c). (18) 

i 

Here the second and third hues of the right-hand side are, respectively, the s-wave and p-wave pairing terms, and 
a vortex located on the center of the system is incorporated into the phase 4>i of the gap functions. To obtain the 
vortex core states for H, we introduce the Bogohubov quasiparticle operator, 

7^ = 5I["t(*)cIt +^'T(^)ciT +^^i(^)cii]- (19) 

i 

The Bogoliubov-de Gennes (BdG) equations is derived from the relation [H, 7^] — E^^ . We solve the BdG equations 
numerically, and calculate the density profile of quasiparticles for low energy states. 

We, first, consider the case with a purely p-wave pairing. In Fig. 7, the density of quasiparticles plotted on the xy- 
plane is shown. There are low energy vortex core states. The lowest energy for the vortex core state is i? 0.000767t. 
For our choice of parameters in this calculation, A'^/Ep 0.005. Thus, the vortex core state has much smaller energy 
than that of the conventional Caroli-de Gennes-Matricon mode. Furthermore, we find edge states with low energies 
as shown in Fig. 7. The existence of low energy edge states is a concomitant of zero energy modes in vortex cores, and 
characterizes topological nature of the state. Therefore, we conclude that the zero energy vortex core states exist. 

The zero energy vortex core states and the low energy edge states survive even when a Zeeman field along the 
z-direction is applied. The calculated results in this case is shown in Fig. 8. 

We, next, consider the case with the admixture of s-wave pairing and p-wave pairing. The numerical results shown 
in Fig. 9 demonstrate that the topological phase with gapless edge modes and zero energy vortex core states is stable 
as long as the gap of p-wave pairing is larger than that of s-wave pairing. As the magnitude of |As| approaches |At|, 
the gap in one of two SO split bands |At| — |As| decreases, and the spectrum gap in the vortex core becomes smaller. 
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FIG. 8: The density of quasiparticles for the purely p wave state plotted on the 81 x 81 xy-plane for low energies. The top 
and bottom panels are, respectively, the plots of |u|(r)|^ + |u|(r)|^ and |?;|(r)|'^ + |t;|(r)|'^. a = t, At = 0.05t, fisHz — 0.04t, 
fi = ~4.0t. E = 0.00123t (left), E = 0.0182t (middle), E = 0.0315t (right). Edge states appear for E = 0. 01824 and 
E = 0.0315*. 



which makes it difficult to clarify the existence of zero energy vortex core states numerically. The numerical approach 
in this section is applicable only when |At| — |As| is not much smaller than |At|. 

To discuss whether the non-Abelian statistics of vortices is possible or not, we need to examine the degeneracy of 
Majorana modes. The numerical analysis presented in this subsection is not suitable for this purpose because of the 
limitation of the system size used in the numerics. To attack this issue, we exploit an analytical method in Sec. IV D. 



B. A comment on the use of the chirality basis for the issues of vortex core states 



Here, we discuss the limitation of the use of the chirality basis for the issues of vortex core states. For concreteness, 
we consider a noncentrosymmetric s-wave superconductor. In the absence of Zeeman fields, in homogeneous bulk 
systems without vortices, the unitary transformation (afc+,afc_)"^ = C/(A;)(cfe-f , Cfc^)^ with 

(20) 

which diagonalizes the asymmetric SO interaction in the spin space transforms the s-wave pairing terms into the 
intra-band pairing terms in the chirality basis. 




5:5:A.e-«('=)aLaL,.+h.c., (21) 



v=± k 



where e^'^*^'^' = —{£ox ± '''^Oy) I \j ^ox + '^Oy odd-parity phase factor. The above expression is similar to the 

pairing term oip-\-ip superconductors. Thus, one may expect that zero energy vortex core modes exist when vortices 
are introduced into the system. However, this is not true. The unitary transformation (20) is derived for spatially 
homogeneous systems, and generally its use is not allowed in the case with vortices. In the case that quasiclassical 
approximation is justified, and the center of mass coordinate of Cooper pairs which characterizes inhomogeneity of 
the system is approximately commutative with the momentum operator k and the unitary transformation U{k), the 
chirality basis can be used for the description of the vortex state, as done in [53-55]. However, to examine whether 
the zero energy mode exists or not in a vortex core, one needs to go beyond the quasiclassical approximation; i.e. the 
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FIG. 9: The density of quasiparticles for the s + p wave state plotted on the 81 x 81 xy-plane for low energies. The top and 
bottom panels are, respectively, the plots of |w|(r)|^ + |u|(r)p and |w|(r)p + |u|(r)p. a = t, At — 0.05t, As = — 0.02t, 
MbHz = 0, fi= -4.0t. E = 0.0023t (left), E = 0.039t (middle), E = 0.0434t (right). Edge states appear for E = 0.0434t. 

unitary transformation for homogeneous systems (20) is not applicable to the issue that we are concerned with here. 
Missing this consideration in [34] leads to an incorrect conclusion that zero energy vortex core states exist even for 
purely s-wave NCS. For the clarification of the existence of zero energy modes in vortex cores, a powerful and reliable 
approach is to exploit the method developed by Tewari, Sarma and Lee which is analogous to the Jackiw-Rebbi's 
index theorem [35, 56]. We consider this issue in the subsequent sections. 

C. Majorana zero energy modes in vortex cores: Index theorem 

In this subsection, we develop an analytical approach for zero energy vortex core states which is a generalization 
of the index theorem for p + ip superconductors obtained in [35] to the case of noncentrosymmetric systems. We 
consider the model for a Rashba superconductor in two dimensions with a vortex located on the center of the system. 
The Hamiltonian is 

Ti- ~ Ti-K + ^so + ^Ps + Ti-pt, 

k.a 

Hpt = -^At I rf'firf're'^«;i,(fi)[gT(r)4+,^4_,^ + <?x(r)4+,^ct,_^J +h.c., (22) 

where the band dispersion in the case with the Zeeman field Hz is ^fco- = ej, — /i — a^L^Hz with tk = k'^ /Ira and /i a 
chemical potential. Tiso is the Rashba SO interaction with the coupling constant a. Tips and Tipt are, respectively, 
the pairing interaction for spin-singiet and spin-triplet channels. The center of mass coordinate and the relative 
coordinate for Cooper pairs are, respectively, represented by R and r. The phase (pR is due to the vortex at i? = 0. 
The suppression of the superconducting order parameter in the vicinity of the vortex is incorporated into the functions 
ha{R) = 1 — e"!-'^!/^" with a = s for the spin-singlet pair and a = t for the spin-triplet pair, gsir), g^ir), and gi{r) 
are the structure functions corresponding to the pairing symmetry. For simplicity, in the following, we consider the 
case with s-wave pairing for the spin-singlet channel and the p-wave pairing for the spin-triplet channel. Then, the 



14 



Fourier transforms for the gap structure functions are gs{k) = 1, = —i{kx — iky) and gi{k) = —i{kx + iky). 

Here the form of the d-vcctor for triplet pairs are determined so as to be consistent with the Rashba SO interaction. 
To discuss the zero energy state in the vortex core, we follow the approach developed in [35], and use the angular 
momentum representation of the electron operators, 

^ oo 
* m= — cx) 

where k = \k\ and 4>k is the azimuthal angle of fe. Cm^ka satisfies the anticommutation relation {cm,ka,cl^ p^,} = 
SmnSaa'S{k — p)- In this representation, the Hamiltonian is rewritten into the following form. 

7^:32 ^'=<^'=L,feaCm,fea, (24) 

/dk 
(^^^[cL.feT'^^+i.fci ~ cJ„,feiCm-i,feT] (25) 

Hps = -iAs / c!A;dp^w„(fc,p)v^[A:cj;_„ ;-^c|„p^ +pc;L„,fc|c|„+i,pJ + h.c, (26) 

■' m 

Wpt = At / dkdp^[{um{k,p)k-Urn+i{k,p)p}ky/kpc^_^ ,^^c^^.p-^+Um,{k,p)k'^ (27) 

Here Um{k,p) is the Fourier transform of l/|fe 

T^:^ =Y.Urn{k,p)e-^-<^--^-\ (28) 

In Eqs.(24),(25), (26) and (27), the m = mode with up spin and the m = 1 mode with down spin are decoupled 
from other modes. The pairing interaction terms for these modes are 

In Ti}pl\ the term with ui{k,p) for the tt pairs vanishes because of the symmetric property u\{k,p) = ui{p, k). Thus, 



■•Ft 
we have 



Hp^-' = At j dfc<ipv^^^[^^o(fc,p)co_fc|cJ p^ + ui(A:,p)cJ ,,^cj j,^] + h.c.. (30) 
1. Purely p-wave pairing 

We, first, consider the purely triplet case, i.e. Ag = 0. Then, the Hamiltonian for the m = mode with up spin 
and the m = 1 mode with down spin is 

/dk 

- j dkdp[Atoik,p)cl ,^^cl p^ - Atic^^cj + h.c], (31) 
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with Ato(fc,p) = At\/fcp(fc^ — p'^)uo{k,p)/2, and Ati(fc,p) = —At^/kp{k'^ — p^)ui{k,p)/2. Here we have antisym- 
metrized the spin-triplet pairing terms. In the Fourier coefficients, 

the most dominant contributions are from (j)k — <i>p ±7r, which implies 'Ui(A;,p) < and wo(fc,p) > 0. By using the 
unitary transformation 

Com = C+ {k)ak+ + K- {k)ak- , (34) 
ciM = K-{k)ak+ + C+{k)ak-, (35) 

with 



" ' (36) 



we rewrite the Hamiltonian (31) as, 

^(01) = JdkJ dp{A4k,p)alal + h.c.}] 

+ j dk j dp[A2+{k,p)al^al_+A2-{k,p)al_al++h.c.], (37) 

where 

A±(fc,p) = ±^^4^r+(fc,p) + ^^^^r.{k,p), (38) 

A2±(fc,p) = Ti^^4^.-(fc,P) - i^^^^s^{k,p), (39) 

withr±(fc,p) = C+{k)C+{p)±C.{k)C-{p), s±{k,p) = C+(A:)C-(p)±C-(fc)C+(p), andei") = e,- ^I-,.y^aWTiiIM^■ 
We concentrate on low energy excitations in the vicinity of the Fermi surface, and write k and p a.s k = kF± + q and 
p = kp± +q' where kp^ and kp^ are, respectively, the Fermi momenta for the two SO split bands, and \q\, \q'\ ^ kp±. 
In Eq.(39), the first term is of the order 0{{q — q')'^), because An)_i(k,p) and s^{k,p) are antisymmetric with respect 
to the exchange of k and p. We can neglect this term compared to other terms up to 0{q — q'). Furthermore, from 
Eqs.(32) and (33), we see that |Aio — Ati|/|Ato + Ati| ~ 0{1/R^) where i?c is the system size, and in the limit 
of i?c oo, the second term of (39) is negligible compared to the first term of (38). As a result, the inter-band 
pairing terms can be neglected, and the Hamiltonian (37) is decoupled into two parts corresponding, respectively, 
to contributions from the two SO spHt bands, z/ = + and — . A±{k,p) = A±{q — q') must be odd in q — q' , since, 
otherwise, the pairing term in Eq.(37) vanishes. Then, the Fourier transforms of A±{q — q') denoted as ^im±{x) are 
odd functions of x. Here m±{x) is real. 

Linearizing the band dispersion around the Fermi momentum, and expressing the operator for the Bogoliubov 
quasiparticle as, 

7^ = j dx[riiual{x) +r]2iyat,{-x)], (40) 

we obtain the Bogoliubov-de-Gennes (BdG) equations from the relation [7{^°^\j^ = 

-iv^azdxr]„{x) + um„{x)ayr]„{x) = Er]^{x), (41) 
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with v = zt corresponding to the two SO split bands and the Fermi velocities v,y. Here /7^(x) = (?7i±, '72±)- The BdG 
equations have a zero energy solution with E = for each bands: When m±{x) satisfies m±{x) = sgii{x)\m±{x)\, the 
zero energy solutions are 

,l{x) = e-^^^'--^^'''^il,l), (42) 

e{x) = e-^^o"'-^'^''j=il,-l). (43) 

For these solutions (42) and (43), the quasiparticle operator (40) satisfies ^yl = 7,y, and thus there are two Majorana 
fermion modes corresponding to the two bands. 

It is noted that these two Majorana modes are stable against the magnetic field along in accordance with the 
existence of two gapless edge modes in this case as discussed in Sec. Ill B. This phase is topologically equivalent to 
the phase I or I' in Table I. 

2. s + p wave pairing 

Now we consider the case with the admixture of the spin singlet pairs and spin triplet pairs. The Hamiltonian for 
the m = mode with up spin and the m = 1 mode with down spin is 



7^(01) 

+ ia 



dk * I 

- dkdpAs{k,p)[cl ,.-^c\ p^+ h.c] 

+ j dkdp[Atoik,p)cl^^clp^ - An{k,p)clk^clp^ + h.c.]. (44) 

Here As{k,p) = Asy^kp[ui{k, p)k + uo{k,p)p]. The application of the unitary transformation (34) and (35) to the 
pairing terms of (44) gives, 

j dkdp[A'^{k,p)al^al+ + A'_{k,p)al_al_ + h.c] 

-i J dkdp[iA'2{k,p)al_^_al_ - iA'2{p, fc)a|._a^_,_ + h.c], (45) 

where 

A'±(fc,p) = A±(fc,p) T Ai'\k,p)s-{k,p) - Ai-\k,p)s+{k,p), (46) 

A',{k,p) = Ai''> {k,p)r+{k,p) + Ai-\k,p)r_{k,p), (47) 

with Ai''\k,p) = (As(fc,p) + As(p, fc))/4 and Ai"^(/c,p) = (As(fc,p) - As(p, A;))/4. It is reasonable to postulate that in 
the s-wave gap Ag, the symmetric part Ag"'' dominates, and Ag"'' can be neglected. We write k = kF± + q, p = kF±+q' 
with q,q' <C fc^i as before. A'_^_{k,p) = ±At{q — q') is odd in g — q' , while A'2{k,p) = As{q — q') is even in g — q' . 
Fourier transforming to the coordinate space, we introduce the quasiparticle operator. 



= J dx[r]i+{x)al_{x) + ri2+{x)a+{~'x) + r]i_{x)ai{x) + ri2-ix)a-{-x)]. (48) 
The Fourier transform of the odd-parity pairing term is 

u J dx imt{x)al{x)al{—x) +h.c, u = ±, (49) 
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while for the even-parity pairing term, 

i/d.m.(.)at +h.c.. (50) 

Here, mi{x) (rns(x)) is the Fourier transform of Ai{q) {As{q)) and odd (even) in x. Then, the BdG equations for 
*^(a;) = {r]i+{x),r]2+{x),r]i_{x),r]2-ix)) are, 

-ivd^r]i+{x) - imt{x)r]2+{x) - ims{x)r]2- (x) = Eiji+ix), 
ivdxr]2+{x) + imt{x)r]i+{x) - ims{x)r]i_{x) = Er]2+{x), 

-ivdxrii-{x) + imt{x)ri2-{x) + ims{x)ri2+{x) = Er]i^{x), 
ivdxT]2-{x) - imt{x)r]i-{x) + ims(x)??i+(x) = Eri2^{x). (51) 

Here, to simphfy the analysis, we have assumed t;+ « f_ = f which is justified for Ep ^ a. From Eq. (51), we find 
two sets of solutions with zero energy eigen value -E = up to normalization factors, 

*f(a;) = (C+(x),C+(x),-C_(x),C_(x)), (52) 



i^l{x) = {-iC-{x),-iC-{x),iC+{x),-iC+{x)), (53) 

with 

C±(x) = e~" -I- g~ " J^^- '^!/[™t(3/)-TOs(y)] ^g^^ 

and x± the solution of mt{x±) ± ms{x±) = 0. We can easily verify that X- = —x+ and C+{x) (C_(x)) is an even 
(odd) function of x. The above solutions (52) and (53) arc normalizablc only when mt(.x) ± ms(.x) > for x > x± 
and mt{x) ± ms{x) < for x < x±. Therefore as long as the gap for the spin-triplet pairs is larger than that for 
the spin-singlet pairs, the zero energy modes exist, which is consistent with the recent result obtained by Lu and Yip 
[38]. It is noted that the quasiparticle operator (48) for the solutions (52) and (53) satisfies 7''' = 7; the quasiparticles 
corresponding to these solutions are Majorana fermions. 

The phase with the two Majorana fermion modes is classified as the phase I or I' in Table I. 



D. Non-Abelian statistics of vortices 



The non-Abelian statistics of vortices is realized when there is only one Majorana mode in a vortex core [8, 9, 17]. 

Thus, it is necessary to eliminate one of two Majorana fermion modes found in Sec. IV C. For this purpose, we 
consider the case that the Fermi level crosses the F point in the Brillouin zone; i.e. £fe=o — /U = (/U = — 4t). Then, 
for ^ 0, a gap ~ /XB-ffz opens in the vicinity of the T point at the Fermi level [34]. 

In the case of purely wave pairing, this implies that u_ in (41) vanishes, and instead the mass term 2(7z\iJ-BHz\'ri-{x) 
is added. In this case, there is no zero energy mode for the v = — band, and there is only one zero energy Majorana 
mode for the quasiparticles with the Fermi momentum kp+ij^ 0) in the ly = + band, which ensures the non-Abelian 
statistics of vortices [34]. This state is topologically equivalent to the phase II in Table I, and also to spinless 
p + ip superconductivity. However, the realization of this state in NCS is more feasible than that of spinless p + ip 
superconductivity, because, for spinless p + ip state, the strong magnetic field associated with full spin polarization 
leads to the fatal orbital depairing effect on superconductivity, while, for the Rashba NCS with /x = —At, a weak 
magnetic field between i?ci and Hc2 applied parallel to the 2;-axis is sufficient to eliminate one of two Majorana 
modes. 

In a similar manner, in the case of s + p-w&ve pairing with |At| > |As|, for ;U = —At, there remains only one zero 
energy mode: The BdG equations for the quasiparticle operator (48) with E = becomes 

-ivdxrii+ix) - imt{x)r]2+{x) - ims{x)ri2-ix) = 0, 
ivdxri2+{x) + imt{x)r)i+{x) - ims{x)r]i-{x) = 0, 

\2fj,BHz\-r]i-{x) +imt{x)rj2-{x) + ims{x)ri2+{x) = 0, 
-\2fj,BHz\m-ix) - imt{x)T]i_{x) + ims{x)T]i+{x) = 0. (55) 

When mt{x) satisfies mt{x) = sgn(a;)|TOt(x)|, (55) has only one normalizable solution, 

t ^ rit \ nf \ t \ "^s(a;)(mt(x) - 2i|^B^^z|) . ^ ms{x){nit{x) - 2i\fj.BH,\) 
m^{x) = C{x), m^ = C{x), m-W= 4(;.Bg.)^+mg(.) ' ^^"^"^ = Ai^BH^r+m^x) ' ^^'^ 
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with 



C{. 



a;) = e " 



4((<BH^)^+">t(!') 



For this solution, the quasiparticle operator is 



7 



= j dxC{a 



4(MB//.)2+m2(a;) 



a+ix) + a+{-x) + 



(i_{x) — a-{—x)) 



(57) 



(58) 



which satisfies 7^ = 7. Thus, there is only one Majorana zero energy mode in a vortex core. Under this situation, the 
non-Abehan statistics of vortices is possible. 



E. Majorana condition 

The Majorana condition of zero energy vortex core states is crucial to the non-Abelian statistics of the vortices, so 
it is better to argue it without any approximation. In this section, we present a general argument on the Majorana 
condition of vortex zero modes for the 2D NCS. 

Let us start with (44) in the Nambu representation. 



(Qi) 1 f dkdp /ft \ / \ 

2 J l^o,fcT''^i,fei''=0''=T.Ci,/ciJ H(fc,p) 



where Hlkjp) is given by 



( ^k'\5k,p iak5k,p h^(k,p) h°^{k,p)\ 

-iakSk,p CkiSk,p -hf-{p,k) h°l{k,p) 

h^^*{p,k) -ft.°^*(fc,p) -^k'\Sk,p iakSk,p 

\h°'*{p,k) h°l*{p,k) -iak5k,p -^kiSk,pJ 



with 



h°Hk,p) = -i{2nfA,{k,p), h°l{k,p) = 2{2nfAto{k,p) h°l{k,p) = -2{2nfAti{k,p). 
Using the following relation. 



/ Co,pT \ 



= r 



'^o,pT 

V <Pi J 



l2x2 
l2x2 



-l,pl / \ -\.,pi 

one can show that l-iik^p) has the so-called particle hole symmetry, 

rH(A;,p)r = -H*(fc,p), 

which also can be checked directly from (60). 

In the momentum space, the BdG equation is given by, 



y" dpn{k,p)u{p) = Eu{k), 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 



where u{p) is a four component vector u{p) = {a{p), (3{p),'y{p),6{p)y . When u{p) is a solution of the BdG equation 
with the energy E, then, using (63), we can show that ru*{p) is a solution with the energy —E, 



j dpn{k,p)Tu*{p) = -Eru*{k). 



(65) 
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Therefore, if tto(p) is a zero mode of the BdG equation, then Tu'^ip) is also a zero mode. This means that if the BdG 
equation has only one independent zero mode Uo(p), then ito(p) and Tuo{p) should not be independent of each other. 
In general, if there are an odd number of independent zero modes, then at least for one solution, Uo{p) and Tuq{p) 
are not independent. 

For simplicity, suppose that the BdG equation has a single zero mode uo{p)- As mentioned above, tio(p) and ruQ{p) 
are not independent, so we have 



uoip) = Tuo{p), 

by choosing a suitable phase of uo{p). In terms of the components of uo{p), this becomes 



uo{p) 



(66) 



ao(p) ^ 






Poip) 




VoiP) 


Xoip) 




a*{p) 


\ Vo{p) J 




Uo(P)/ 



(67) 



To obtain the annihilation operator 7 for the zero mode, we perform the mode expansion for {copi,Cipi, cj^i, c\p^), 



(68) 



^ co,pT ^ 




^ aoip) \ 




= 7 


Mp) 




Xo{p) 


V 4m J 




\ Vo{p) J 



where we have omitted the terms including non-zero modes. In this equation, we have two dependent relations, 

Co.pT = 7Q!o(p) + • • • , 4,pT = 7Xo(p) + • • • , 
ci,pi = 7/?o(p) + • • • , c\p^ = 7%(p) + • • • • 
Thus 7 must satisfy the Majorana condition, 

7 = 7'l'. 

From the normalization condition for the zero mode, 7 can be written as 

7 = y dp{a*o{p),(3*o{p),x*o{pUo{p)) 



(69) 

(70) 



Cl.pi 

Cq.pT 



(71) 



Thus the commutation relation of 770 can be calculated as 

{7,7'^} = j *M"o(^)co,fcT +/3o(^)ci,fei +Xo(^)co,feT +^o(^)4,fci. 



/ 



"o(p)4,pT + i^o{p)c\^pi + Xo(p)co,pT + Vo{p)ci,pi} 
dkdp [a*o{k)ao{p) + l3*o{k)l3o{p) + Xo{k)xo{p) + Vo{k)mip)] Sk, 



= 1 



(72) 



In a similar manner, we can show that 7 anti-commutes with annihilation and creation operators for non-zero modes. 

In conclusion, the existence of only one zero energy mode in a vortex core is the necessary and sufficient condition 
for the existence of a single Majorana mode which leads to the non-Abelian statistics. 



V. SUMMARY 



We have explored topological phases of NCS characterized by the existence of gapless edge states and Majorana 
fermion modes in vortex cores, mainly focusing on the 2D Rashba superconductors with the admixture of s-wave 
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pairing and p-wave pairing. It has been clarified that when the p-wave gap is larger than the s-wave gap, topological 
states are realized. We have completed the topological classification of these states by examining the TKNN number, 
the Zi index, and the winding number associated with specific symmetry points in the Brillouin zone. It has been 
also found that for the Rashba superconductors, the gapless edge states that ensure the quantum spin Hall effect 
protected from disorder are stable against a weak magnetic field applied perpendicular to the propagating direction of 
the edge states, despite broken time-reversal symmetry which fiaws the Zi characterization of the topological phase. 
The stability of the edge states is guaranteed by a specific accidental symmetry of the Rashba model. We have also 
proposed a simple scheme for the realization of the non-Abelian statistics of vortices in topological phases of NCS 
under an applied magnetic field. 

Experimental verification of these findings are particularly of interest. The topologically-protected gapless edge 
states play important roles for transport properties. In the superconducting state, current flows carried by the edge 
states can be detected by the measurement of a spin Hall current or thermal transport measurements at sufficiently 
low temperatures where bulk quasiparticles are well suppressed. Also, the accidental topological phase in the case 
with a magnetic field can be detected by observing the dependence of the transport current on the direction of the 
magnetic field, or the splitting of the zero-energy bias peak of a tunneling conductance due to the tilt of the magnetic 
field. More elaborate argument on these experimental implications will be addressed in the near future. 

The experimental realization of the non-Abelian statistics of vortices is most challenging, though the scheme pro- 
posed in this paper is, in principle, feasible. The s -hp- wave NCS with the Fermi level tuned to be /U = — 4t discussed 
in Sec. IV D need not to be bulk systems. A proximity- induced superconductor realized in the vicinity of the interface 
between a p-wavc superconductor and a semiconductor with the Rashba SO interaction may be a good candidate for 
the realization of this phenomenon. For the experimental detection of the non-Abelian statistics, the two-point-contact 
interferometer considered in Refs.[5, 10, 34, 57] may be employed. 
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